Stability and dynamics of free magnetic polarons 



M. J. Calderon^'^^ l Brey^ and P. B. Littlewood^ 
^ Instituto de Ciencia de Materiales ( CSIC). Cantoblanco, 28049 Madrid. Spain. 
^ Departamento de Fisica Teorica de la Materia Condensada, Universidad Autonoma de Madrid, 28049 Madrid. Spam. 
^ Cavendish Laboratory, Cambridge University, Madingley Road, Cambridge CBS OHE, UK 

(February 1, 2008) 

The stability and dynamics of a free magnetic polaron are studied by Monte Carlo simulation of 
a classical two-dimensional Heisenberg model coupled to a single electron. We compare our results 
to the earlier mean-field analysis of the stability of the polaron, finding qualitative similarity but 
quantitative differences. The dynamical simulations give estimates of the temperature dependence 
of the polaron diffusion, as well as a crossover to a tunnelling regime. 
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Recently, the interest on free magnetic polarons (FMP) 
has been renewed due to its relation with colossal mag- 
netoresistance [0 in manganese pyrochlores and in 
double-exchange models of magnetism in the manganite 
perovskites Q-^. 

Free magnetic polarons are magnetically self-trapped 
carriers in contrast to the more common bound magnetic 
polarons which become trapped by an impurity. These 
were extensively studied in diluted magnetic semiconduc- 
tors [0 and rare earth chalcogenides ||]. Golnik et al. 

found experimental evidence of the existence of not 
only bound but free magnetic polarons in Cdi-xMuxTe 
and Pbi-xMuxTe. Their existence leads to an activated 
behavior of the resistivity above the Curie temperature 
that is found in materials with negative magnetoresis- 
tance, as well as a temperature- dependent spin-splitting 
seen in magneto-optical experiments pO[ |. Theoretical 
models of FMP have been developed within a mean- 
field approach |pT|- ^3| and generalized to a fluctuation- 
dominated regime |14| . In most of these systems, the 
underlying (super)-exchange interaction between the lo- 
calized spins is antiferromagnetic in nature; we shall how- 
ever be concerned with the ferromagnetic case. 

We study the model of a single electron interacting 
with a spin background that itself is ordering ferromag- 
netically We consider the Hamiltonian: 

H =-t C^Cjcr 

-Jj2{il~ a)(5f • + Sf ■ S]) + St ■ S^} 
-J'Y.3,-S, (1) 
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where Si refer to the spin of the magnetic ions in the 
system. creates an electron with spin a on the site i, 
~ '^ia'^af3Cip IS the couduction spin operator and {i,j) 
denotes sum over the nearest-neighbors pairs. We have 
added to the Heisenberg term a small Ising anisotropy 
a = 0.1 to improve convergence at low temperatures 
and by enforcing a non-zero transition temperature Tc 
[p^ . J' is the coupling between a localized spin and 



a conduction electron. The qualitative behavior is well 
understood from previous mean-field analyses |jl|,|ll|,|l3| . 
Below a temperature Tp a ferromagnetic polaron forms 
by self-trapping in a ferro-magnetically aligned cluster of 
spins. As the temperature is lowered toward the Curie 
temperature Tc the polaron grows in size and becomes 
more stable, because the small-q magnetic susceptibility 
is growing. Near and below Tc the polaron will again 
become unstable because of the ease of motion in the 
background ferro-magnetic spin alignment. Notice that 
this is quite different from the case of anii-ferromagnetic 
coupling of spins, where the polaron may remain stable 
well below the magnetic ordering temperature. 

There are several deficiencies of the mean field treat- 
ment. The most pronounced is a continuum treatment 
of the spin background where fluctuations are neglected. 
This approximation is such that the paramagnetic state 
leads to a vanishing exchange coupling, so that the elec- 
tron is bound in a potential of depth J'S, with S the 
average magnetization inside the polaron. As J' /t — > oo, 
the potential well becomes arbitrarily deep. This is un- 
doubtedly a severe overestimate. 

In the paramagnet there will always be low energy 
states localized in the band tail with energies 0{t) 
above the ferromagnetic ground state, even in the strong 
coupling limit. Such low energy states are produced by 
random fluctuations of a few neighboring spins into near- 
alignment. But now one must distinguish between a self- 
trapped polaron and a localized band-tail state, if indeed 
such a distinction is appropriate. 

In this paper we address the topic by a dynamical sim- 
ulation of the Hamiltonian of Eq.(l). We show that po- 
larons may be distinguished (when they exist) by a spec- 
troscopic gap to band-like states, and that they move 
diffusively. As temperature is raised, the polaron level 
moves toward the band edge, and begins to resonate with 
states in the band tail, leading to a crossover to hopping 
conductivity. 

We perform a classical Monte Carlo simulation (MC) 
in two dimensions on a square lattice of localized spins Si 
which are treated as classical rotors characterized by the 
angles 9i and 0^. We use periodic boundary conditions 
on a two dimensional lattice of size up to 30 x 30. 
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We place a single electron in the system in the lowest 
energy eigenstate of the Hamiltonian consisting of the 
first and third terms of Eq (1), using the instantaneous 
spin configuration for the classical spins Si . The resulting 
wave function leads to a local magnetic field proportional 
to \^{x,y)\'^, used in the next step of the MC spin simu- 
lation. 

The standard Metropolis algorithm is used. Randomly 
chosen sites suffer a random change of spin orientation. 
Changes are allowed if the increment in energy AE is 
such that the quantity exp{—AE/KT) is smaller than a 
random number between and 1. 4000 reorientations per 
spin were made for an initial equilibration and 3000 to 
calculate averages after each diagonalization. Each diag- 
onalization defines our time step. Changing the number 
of spin reorientations between each diagonalization led 
to no significant change in either the magnetization or 
binding energy. All the quantities are given in units of 
J, the Heisenberg parameter. The hopping parameter t 
is fixed to 100 (estimated with the mean field relation 
Tc ~ zJS^ and the values for the parameters expected 
for the pyrochlores j^]) as we are interested in the behav- 
ior versus J' /t and the temperature T. Tc — 1.8 in these 
units. 

This new approach allows us to calculate in a self- 
consistent way the wave function and the magnetic polar- 
ization over a large range of temperature. In particular, 
we can explore the region around Tc where the mean-field 
treatment fails. 

In Fig. |l| we plot \ilj{x,y)\'^ and the averaged local 
magnetization close to Tc- Visually, the existence of a 
magnetic polaron is clear, and there is substantial align- 
ment of the moments in the vicinity of the carrier. Note 
that far from the influence of the wave-function the av- 
erage magnetization is close to 0, so a large part of the 
spin configurations are explored, while for the spins close 
to the center of the wave-function there are few accepted 
spin fiips. 




FIG. 1. |'!/'(a;, and the averaged magnetization 

mi = jSj for J'/t = 5 and T = l.lTc are plotted. 

Pictorial evidence is purely qualitative, and does not 



allow one to extract reliable estimates for the polaron size 
or binding energy, especially at higher temperatures and 
lower J' /t, when the polarons are smaller and fiuctuating 
in time. More reliable evidence comes from the time- 
averaged electronic density of states (DOS), and of the 
excitation spectrum shown in Fig. |^. In the density of 
states (inset) a sharp low energy feature is pulled from the 
bottom of the band (only the lowest 1% of the spectrum 
is shown) that contains exactly one state. This is the 
bound polaron level. The level width comes from thermal 
fluctuations in the energy of the bound state, and the 
stability of the bound polaron is seen more clearly in 
the excitation spectrum (main figure) that demonstrates 
a clear gap corresponding to the electronic part of the 
binding energy Ep of the polaron. 
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FIG. 2. The density of states (inset) and excitation spec- 
trum at low energies for J' /t = 5, T — l.lTc and \S\ = 3/2 is 
plotted. 

The continuum of excited states can be characterized 
as the band tail formed by the fluctuating paramagnetic 
background; the lowest energy states are produced by 
rare fluctuations of nearby spins into near- alignment. 
Consequently, the "gap" in the excitation spectrum is 
soft, and indeed statistically very rare states may occur 
at energies below the bound state of the polaron. We will 
discuss this below. 

We estimate the electronic binding energy Ep by the 
configurationally averaged gap A = £'i — to the lowest 
excited state in our simulations (we have checked that the 
separation between excited states scales as so that 
these are true continuum states). In Fig. || we show 
the dependence of A and the absolute value of the local 
magnetization M (weighted with the wave-function) on 
T and J' /t. M is defined as 

M^{\J2SI\) 
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where §1 = ^. 

We are not taking into account thermal excitations of 
the quasiparticle so our results are valid only when A is 
bigger than T. This condition is fulfilled for all the values 
of A shown in Fig. As expected from previous analy- 
ses, the polaron binding energy increases as temperature 
is lowered from high temperatures, as the thermal spin 
fluctuations are reduced. For large J' /t we find a new be- 
havior on A not found within mean-field theory, namely, 
that it has a maximum at some temperature above Tc- 
The existence of a maximum can be understood in terms 
of the correlation length ^. This quantity increases as 
we decrease the temperature above T^. For very small 
^ (large T) is very difficult to have a FM cluster for the 
spin-polaron to sit in and for large ^ the electron would 
rather spread out. This will lead to an intermediate opti- 
mum 1^ for the existence of the polaron that would happen 
close to Tc but not necessarily at Tc- 
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FIG. 3. We plot here the dependence on J' /t and tempera- 
ture of the local magnetization in the lattice produced by the 
polaron (a), its binding energy (b) and the mobility fi = j. 
(c) for |S| = 3/2 . The different curves correspond to differ- 
ent J' /t such that its value is 1 for open circles, 1.5 for closed 
circles, 2 for open squares, 3 for closed squares, 4 for open 
diamonds and 5 for closed diamonds. 

The size of the polaron may be estimated from the sep- 
aration between the first eigenvalue Eo and the bottom of 
the band of the uniform ferromagnet. The bottom of the 
band in this case is given by — | J' — 4t and the separation 
should go roughly as being Lp the size of the polaron 

if we assume saturation in the local magnetization. The 
general trend is that the size decreases as T (above T^) 
or J' increases. From Fig. ^ we can also deduce that the 
'window' above Tc where the spin-polaron is stable in- 



creases with J'/t. These two results are consistent with 
previous mean-field calculations . 

Although qualitative comparison is satisfactory there 
are large quantitative differences that point to a great de- 
creasing in the stability of the spin-polarons when fluctu- 
ations are taken into account. To be precise we compare 
the binding energy at T — Tc- From mean- field calcu- 
lations on ref. j^] the maximum possible value for A is 
^ J' S but it is not reached due to the kinetic energy that 
is lost with the formation of a polaron. In the present 
work A/ J ~ 1 while J'S/J ^ 100. So binding energies 
are reduced by two orders of magnitude compared with 
the mean-field results because the loss of kinetic energy 
is not well taken into account in the latter. 

The study of the stability conditions for a free magnetic 
polaron is interesting by itself but the MC simulation also 
opens us the possibility of learning about its dynamics in 
a spin-fluctuating landscape. In Fig.^ the probability of 
moving a distance r (defined as the change in the ex- 
pectation value of the electron position) for different MC 
times is shown. For time t = 1 one observes dominant 
short distance motion with occasional rare hops over long 
distances. For longer times, the peak of the distribution 
moves out approximately with 's/I as expected. This is 
the expected behavior from a diffusing object. The long- 
distance hops occur when unoccupied band tail states 
(which may be localized anywhere in the system) tem- 
porarily drop below the bound polaron level. In our algo- 
rithm - which automatically populates the lowest energy 
level - the electron moves to occupy this new state and 
restabilises the polaron there. These rare events even- 
tually dominate the long-time behavior in our simula- 
tions. Of course, very long range hops are unphysical 
because the tunnelling probability will be exponentially 
small with distance, and the band-tail states survive in 
one place for only a short time. Hops to band-tail states 
will then be limited to some finite range. As tempera- 
ture is raised, and Ep is reduced, hops to band tail states 
become more frequent; we cross over to a regime of "pas- 
sive advection" of the wavefunction in the fiuctuating 
spin background [0. 

Our results are fitted to a gaussian in two dimensions 
plus a constant (to approximately take account of hops 
to band-tail states). The gaussian dominates for the pa- 
rameters of interest when a spin-polaron is well formed. 
The distribution scales with ^/t as expected for diffu- 
sive motion. Hence we calculate the diffusion constant 
as D = P(r, t = l)r^ and the mobility (/i — ^) of the 
spin polaron for different couplings and temperatures (see 
Fig. ^). The mobility decreases with temperature, and 
also with J'/t. The latter is reasonable, because larger 
polarons should diffuse more slowly. The temperature- 
dependence is more surprising, and arises because D it- 
self is weakly T-dependent. Although the polaron size 
is decreasing with temperature (tending to increase Z?), 
this is counterbalanced by a reduced probability of favor- 
able FM spin configurations near its boundary as T/Tc 



3 



is increased. 
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FIG. 4. The probability of moving a distance r for different 
times, t = 1, 10 and 100, taken from a single run, is shown. 
Solid lines show fits to a 2D Gaussian, plus an offset. This 
background is due to the rare appearance of FM clusters away 
from the polaron location, where the polaron hops a large 
distance. The curves for f = 10 and t = 100 are the result 
of iterations of the t = 1 curve; thus the rare long-distance 
hops eventually dominate the distribution, as can be seen in 
the trace at the longest times. 

The Heisenberg term has been considered ferromag- 
netic to compare with the pyrochlores. The change to 
antiferromagnetic coupling is straightforward and in fact 
the more common case in manganite perovskites , rare 
earth chalcogenides ||], or magnetic semiconductors 0. 
We find in the case of an antiferromagnetic background 
the stabiHty of a free magnetic polaron is enhanced, and 
will report these results elsewhere. These results are con- 
sistent with ref. where a pseudogap in the DOS is 
associated with phase separation, that is the large scale 
effect corresponding to spin-polarons. 

In conclusion, our dynamical simulations have revealed 
a picture of the FMP in a ferromagnet above Tc which 
is considerably more complex than given by the mean 
field pictures. Provided the exchange coupling is large 
enough, FMP's are stable above Tc, but considerably 
more weakly bound than found by mean field calcula- 
tions. This by itself raises some doubts about the inter- 
pretation given earlier for the Mn pyrochlores, because 
we require an exchange coupling comparable to the band- 
width for a well-formed polaron with nearly saturated 
magnetization, whereas in the Mn pyrochlores this cou- 
pling is expected to be not large Q. We find that the 
motion of the polaron is diffusive, but as temperature is 
raised the electron fluctuates out of the self-trapped con- 
figuration into band-tail states formed by opportunistic 
fluctuations of the moments. 
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